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1. Introduction 

Cell membranes play crucial role in living movements. They consist of lipids, proteins 
and carbohydrates etc. There are many simplified models for cell membranes in 
history [1]. Among them, the widely accepted one is the fiuid mosaic model proposed 
by Singer and Nicolson in 1972 [2]. In this model, a cell membrane is considered as a 
lipid bilayer where lipid molecules can move freely in the membrane surface like fluid, 
while proteins are embedded in the lipid bilayer. This model suggests that the shape 
of the cell membrane is determined by its lipid bilayer. Usually, the thickness of lipid 
bilayer is about 4 nanometers which is much less than the scale of the cell (about several 
micrometers). Therefore, wc can use a geometrical surface to describe the lipid bilayer. 
In 1973, Helfrich [3] proposed the curvature energy per unit area of the bilayer 



spontaneous curvatures of the membrane surface, respectively. We can safely ignore 
the thermodynamic fluctuation of the curved bilayer at the room temperature because 
of kc ~ 10~^^J ^ ksT [4, 5], where is the Boltzmann factor and T the room 
temperature. Based on Helfrich's curvature energy, the free energy of the closed bilayer 
under the osmotic pressure p (the outer pressure minus the inner one) is written as 



where dA is the area element, // the surface tension of the bilayer, and V the volume 
enclosed within the lipid bilayer. Starting with above free energy, many researchers 
studied the shapes of bilaycrs [6,7]. Especially, by taking the flrst order variation of the 
free energy, Ou-Yang et al. derived an equation to describe the equilibrium shape of the 
bilayer [8]: 



They also obtained that the threshold pressure for instability of spherical bilayer was 
Pc ~ kc/R^, where R being the radius of spherical bilayer. 

Recently, opening-up process of lipid bilayers by talin was observed by Saitoh 
et al. [9,10], which arose the interest of studying the shape equation and boundary 
conditions of lipid bilayers with free exposed edges. Capovilla et al. flrst gave the shape 
equation and boundary conditions [11] of open lipid bilayers. They also discussed the 
mechanical meaning of these equations [11,12]. In recent paper, we also derived the 
shape equation and boundary conditions in different way — using exterior differential 
forms to deal with the variational problems on curved surfaces [13]. It is necessary to 
further develop this method because we have seen that it is much more concise than the 
tensor method in recent book [6] by one of the authors. 

In fact, the structures of cell membranes are far more complex than the fluid 
mosaic model. The cross-linking structures exist in cell membranes where fllaments 
of membrane skeleton link to proteins mosaicked in lipid bilayers [14]. It is worth 





(2) 



p - 2nH + kc{2H + co){2H^ - cqH - 2K) + kcW\2H) = 0. 



(3) 
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discussing whether the cross-hnking structures have effect on the shapes and stabihties 
of cell membranes. 

In the following contents, both lipid bilayers and cell membranes are called bio- 
membranes. We will fully develop our geometric method to study the shapes and 
stabilities of bio-membranes. Our method might not new for mathematicians who are 
familiar with the work by Griffiths and Bryant et al. [15, 16]. Our method focuses on 
the application aspect, but the work by Griffiths and Bryant et al. emphasizes on the 
geometric meaning. Otherwise, we notice a nice review paper by Kamien [17], where 
he give an introduction to the classic differential geometry in soft materials. Here 
we will show that exterior differential forms not mentioned by Kamien might also be 
useful in the study of bio-membranes. This paper is organized as follows: In Sec|21 we 
briefly introduce the basic concepts in differential geometry and the variational theory 
of surfaces. In SecEl we deal with variational problems on a closed surface, and derive 
the shape equation of closed lipid bilayers, and then discuss the mechanical stabilities of 
spherical bilayers. In SecEJ we deal with variational problems on an open surface, and 
then derive the shape equation and boundary conditions of open lipid bilayers as well 
as two-component lipid bilayers. In SecEl we derive the free energy of the cross-linking 
structure by analogy with the theory of rubber elasticity, and regard the free energy 
of the cell membrane as the sum of the free energy of the closed bilayer and that of 
cross-linking structure. The shape equation, in-plane stain equations, and mechanical 
stabilities of cell membranes are discussed by taking the first and second order variations 
of the total free energy. In SecIHl we summarize the new results obtained in this paper. 



2. Mathematical preliminaries 

Here we assume that the readers are familiar with the basic concepts in differential 
geometry, such as manifold, differential form and Stokes theorem (see also Appendix A). 



2.1. Surfaces in three-dimensional Euclidean space, moving frame method 

At every point P of a smooth and orientable surface M in three-dimensional Euclidean 
space , as shown in Fig^ we can construct an orthogonal system ei , e2 , such that 
es is the normal of the surface and e^ ■ e^ = 6ij, {i,j = 1, 2, 3). We call {P; ei, e2, e^} a 
moving frame. For the point in curve C, we let ei be its tangent vector and e2 point to 
the inner point of M. The difference between two frames at point P and P' (which is 
very close to P) is denoted by 

dr = lim PP' = cuiei + ujoeo, (4) 
p_p, ^ ^ 

dei = ujijej (i = 1,2,3), (5) 

where Ui, uj^ and cUjj (i, j = 1,2,3) are 1-forms. 

It is easy to obtain lOij = —uji from e^ ■ ej = 6ij. Additionally, using ddr = and 
ddei = 0, we obtain the structure equations of the surface: 

dui = iUi2 A 102] (6) 




duJ2 = UJ21 f\ UJi] (7) 
UOi A CJi3 + to'2 A to'23 = 0; (8) 

duij = LUik A Wfcj (z, j = 1, 2, 3). (9) 
If we considering the Cartan Lemma, (jH)) suggests 

ciJi3 = acui + hu2 and C1J23 = fecc'i + cuj2. (10) 

Thus we have [18]: 

The area element: dA = oji /\uj2, (11) 

The first fundamental form: / = rfr ■ rfr = cu^ + co'2, 

The second fundamental form: // = —dr ■ de^ = aLuf + 2huJiUJ2 + cu^, 

The third fundamental form: /// = de-^ ■ de^ = ul-^ + 1^32, 

Mean curvature: H = {a + c)/2, (12) 

Gaussian curvature: K = ac — h^. (13) 

2.2. Hodge star * and Gauss mapping 

2.2.1. Hodge star * Here we just show the basic properties of Hodge star * on 

surface M. So as to its general definition, please see Ref. [19]. 

li g,h are functions defined on 2D smooth surface M, then the following formulas 
are valid: 

* / = M A UJ2; (14) 
*uji = 002, = —^i', (15) 
d*df = VV^i A UJ2, (16) 
where is Laplace-Beltrami operator. 
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It is easy to obtain the second Green identity 

/ ifd*dh-hd*df)= [ if*dh-h*df) (17) 
Jm JdM 

through Stokes theorem and the integration by parts. It follows that 

fd*dh= I hd*df, (18) 
M Jm 

if M is a closed surface. 

2.2.2. Gauss mapping The Gauss mapping ^ : M — > S*^ is defined as ^(r) = 

63 (r), where S'^ is a unit sphere. It induces a linear mapping Q* : ^ such that: 

(i) g^uji = ujis, g*uj2 = UJ23; 

(ii) if df = fiuji + /2u;2, then g*df = fiQ*uJi + f2Gy2- 

Thus we can define a new differential operator d = Q*d. Obviously, if df = 

fiUi + f2UJ2, then df = fiUi^ + f2^^23- If define a new operator * such that *u!i3 = UJ23 
and *u!23 = —UJ13, we have 

Lemma 2.1 Jj^^{fd*dh — hd*df) = jQj^j{f*dh — h*df) for the smooth functions f and 
h on M . 

Proof: Using the integration by parts and Stokes theorem, we obtain 

fd*dh= [ f*dh- I dfA*dh, (19) 

M JdM Jm 



hdidf = / hidf - / dhAidf. (20) 

M JdM Jm 

Otherwise, if let df = fiUi + f2UJ2 and dh = hiUi + h2U!2, we can prove 
df A *dh = dh A *df = [0/2/^2 + cfihi — &(/i/i2 + /2^i)]^^i A 002 through a few steps 
of calculations. Therefore, We can arrive at Lemma f2. II bv (fT^ minus (1201) -H 

It follows that 



fdidh= / hdidf (21) 
M Jm 

for a closed surface. 

Because d*df is a 2-form, we can define an operator V ■ V such that d*df = 

V ■ V/cJi A UJ2. 

2. 3. Variational theory of surface 

If let M undergoes an infinitesimal deformation such that every point r of M has a 
displacement 6r, we obtain a new surface M' = {r'|r' = r + 6r}. 6r is called the 
variation of surface M and expressed as 

5r = 5ir + + 83^, (22) 
5.r = ae, (2 = 1,2,3), (23) 

where the repeated subindexes do not represent Einstein summation. 
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Definition 2.1 If f is a generalized function of r {including scalar function, vector 
function, and r-form dependent on point r), define 

5^^f = (g!)/:('^)[/(r + 6,r) - /(r)] {t = 1, 2, 3; g = 1, 2, 3, ■ ■ ■), (24) 

and the q-order variation of f 

5(')/ = (g!K^^H/(r + 5r)-/(r)] (g = 1, 2, 3, ■ ■ ■), (25) 

where £'•''•'[■ ■ ■] represents the terms of ^If^ni^nf' in Taylor series of [■ ■ ■] with qi + q2 + 
qs = q and qi,q2,q3 being non-negative integers. 
It is easy to prove that: 

(i) 6^'^^ and 6'^'^^ {i = 1, 2, 3; g = 1, 2, ■ ■ ■) are linear operators; 

(ii) 5[^\ 5!^\ and ^^^^ are commutative witli eacli otlier; 

(iii) (5f = (5f ^^J^^ and 5(9+^) = thus we can safely replace 6^^\ 
and 5*^^^ by 5i, 5^, 5, and 5' (g = 2, 3, ■ ■ ■), respectively; 

(iv) For functions / and g, Si[f{r) o g{r)] = Sif{r) o g{r) + /(r) o 6ig{r), where o 
represents the ordinary production, vector production or exterior production; 

(v) 5J[g{r)] = {df/dg)S,g; 

(vi) 5" = {6i + 62 + (53)^ e.g. 52 = 61 + 6^ + + 26162 + 26263 + 26163. 

Due to the deformation of M, the vectors ei, 62, 63 are also changed. Denote their 
changes 

6iei = ^Hjej, (26) 

Obviously, ■ ej = 6ij implies Quj = —Vtiji. Because 5i, 62 and ^3 are linear 
mappings from M to M', they are commutative with exterior differential operator d [18]. 



Therefore, using d6iY = 6idv and d6iej = 6idej, we arrive at 








(27) 


61UJ2 = fiico'12 — uJiQiu, 




(28) 


^113 = a^i, ^^123 = ^^1; 




(29) 


62^1 = ^2^21 ^ '-<^2^221; 




(30) 


62LO2 = dVL2 - ^^1^212, 




(31) 


^213 = ^^2; ^223 = cVt2] 




(32) 


63UJ1 = ^230731 — 072^2321, 




(33) 


63UJ2 = ^3'^32 — ^^1^312; 




(34) 


d^3 = ^313^^1 + ^323^^2; 




(35) 


^l^ij = d^lij + ^lik^kj " ^ik^lkj- 




(36) 


Above equations (j2ZI) ~ (jHHjl are the fundamental equations in our paper and have 


not existed in previous mathematical literature 


[15] and [16]. 


Otherwise, it is easy to 


deduce that 6idf = d6if {i = 1, 2, 3) for function 


/• 
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3. Closed lipid bilayers 

In this section, we will discuss the equilibrium shapes and mechanical stabilities of closed 
lipid bilayers. We just consider the closed surface in this section. 

3.1. First order variational problems on a closed surface 

In this subsection, we will discuss the first order variation of the functional 

J^= [ S{2H[T],K[r])dA + p [ dV, (37) 

J M Jv 

where H and K are mean and gaussian curvatures at point r in surface M. p is a 
constant and V be the volume enclosed within the surface. 

According to the variational theory of surface in SecEl we have 5JF = + S2J^ + 
S^J^. Therefore, the next tasks are to calculate SiJ-", 62^-" and S^J-", respectively. 

3.1.1. Calculation of Here, we will briefly prove 4 Lemmas and a theorem. 

Above all, denote 



/ £i2H[r],K[r])dA. (38) 
Lemma 3.1 ^^dA = -{2H)VL-idA. 

Proof: d'idA = 53{uJi A UJ2) = 6-iUJi Auj2 + uji A 53UJ2. Considering (jHSl) and 

we arrive at this Lemma. % 



Lemma 3.2 63{2H)dA = 2{2H^ - K)n3dA + d* d^a. 

Proof: 63{2H)dA = S^auJi A 0^2 + SsOJi A UJ2- Let 63 acts on (fTUj) . we have 

53UJ13 = S^auji + aS^uJi + (53&co'2 + b5-iUJ2, 

If considering (fT^ . (fT!^ . fTHjl . we arrive at this Lemma. ^ 

Lemma 3.3 S-^KdA = 2KHVt3dA + didVt^. 

Proof: Theorem Egregium (see [Appendix C\ implies that S^KdA = —Ssduu— KS^dA = 
—dS^Uu — KS^dA. We will arrive at this Lemma from ()36|) and Lemma 13.11 as well as 
the discussions in Sec J2.2.2l ^ 

Theorem 3.1 ^s^e = JmH^^ + - 2^)a^ + (V ■ V + 2KH)§^ - 2H£]Vt^dA. 
Proof: Above all, we have 

53^e= / 5^8dA+ I 85^A 

JM JM 

= I lJ^^3(2i^)rfA+ / ^5,KdA+ [ £5,dA. 
Jm d{2H) Jm dK 
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+ 



M 



M 



- 2K)J^ + 2ra|^ - 2H£ 



d{2H) ' dK 
a* d\l^ + -——a*a\l3 



d{2H) 



dK 



(39) 



For the closed surface M, we arrive at this theorem by considering (|TH)l and ()2H1 . % 
Lemma 3.4 63 JydV = Vt^dA. 

Proof: Because M is a closed surface in E^, Stokes theorem (see [Appendix A| ) implies 
jy 'idV = Jy'V ■ rdV = Jgy r ■ ndA, thus 



dV = ^ 63[r ■ e3(u;i At^a)] 
3 Jm 



(40) 



From 



(41) 
(42) 
(43) 



= 1 [ l^sr ■ 63(0;! A cua) + r ■ 5363 (lui A c^a) + r ■ 6353 (tUi A 1^2)]- 
3 Jm 

^(jnSI), we obtain 

d^r ■ e3(u;i A 002) = ^3^1 A UJ2, 
r ■ es53{uJi A tua) = r ■ e-i{-2H)Q,-iUJi A cua, 
r ■ 5363(0;! A tua) = - dils A (-r ■ eaCJi + r ■ eiCJa). 
By using the integration by parts and Stokes theorem, we have 

- / dQ^ A (-r ■ eaCJi + r ■ eitua) = / Q^d^-r ■ eacui + r ■ eiCJa) 
Jm Jm 

= I Qsl^ + r ■ e3{2H)]uJi A U2 
Jm 

Therefore, we will arrive at ^3 JydV = j\,jQsdA by using (gOll^dlll)- If 

3.1.2. Calculation of 5iJ^ and 82^ 

Theorem 3.2 5iJ^ = and ^aJ" = 0. 
Proof: We obtain 



(44) 



db AuJi + 2bduJi = (a — c)duj2 — do A c^a- 

from © and (Uni). 

Using (I2I|)~(I2ni), dSni), and (gni), we arrive at 

6i{uJi A UJ2) = (i(fiiCJa), 
6i{2H)ui AuJ2 = d{2H) A tua^i 

through a few calculations. 

By analogy with the proof of Lemma 13. 3| we can prove that 

diKoji Auj2 = dK A fiiCUa- 



(45) 



(46) 
(47) 



(4^ 
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M 



(49) 



M 



Similarly, we can obtain 



d{SuJiQ2)- 



(50) 



M 



Otherwise, it is not hard to obtain 61 jydV = j\,jd{v ■ e^uj2VLi) and 82 JydV 
Therefore, 61J-' = 62J-' = because M is a closed surface. % 



3.1.3. Euler- Lagrange equation 
6J^ = 



Till now, we can obtain 



M 



(V2 + - 2K)^^p— + (V ■ V + 2KH)^ -2HS+p 



d{2H) 



dK 



n^dA.{^l) 



Thus the Euler-Lagrange equation corresponding to the functional is: 



d 



d(2H) 

The similar equation is first found in Ref. [20] 



+ ( V ■ V + 2KH^ ^^'^^ 



£{2H, K)+p = 0.(52) 



3.2. Second order variation 



In this subsection, we discuss the second order variation of functional (j?fj) . This problem 
was also studied by Capovilla and Guven in recent paper [21]. Because = ^2^^ = 
for closed surface M, we have 881^ = 882^^ = 0, and = 88sJ-' = 8lJ-'. 
Form ()39|) and Lemma f3. 41 we obtain 



8^ J" = 8s 



I M 

+ 5z 



BF r)F 
{4H^ - 2K)-j-— + {2KH)^ - 2HF + p 



d£ 
M d{2H) 



d{2H) ' 'dK 

f OF ^~ 
d * dQs + 83 — — d*dQ3 



QsdA 



M 



dK 



8, 



M 



BF BF 
{4H^ - 2K)-j-— + {2KH)-^ - 2HF + p 



+ 



M 



B{2H) ' BK 

BF BF 
{4H^ - 2K)^^^ + C^KH)— -2HF+P 



+ / h 

'm 



BF 



B{2H) 
d * d^l-i + 



fl^dA 



BF 



B{2H) 



_B{2H) 

f BF\ ^~ BF 
+ / ( T777 I d*dVLs + ——8s{d*dVLs). 



8-i{d * dVt-i 



,M \dK J " ' BK"'' ^^^^ 

Please notice that ^3 can freely come into and out of the expressions acted by the 
operator ^3. 
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Lemma 3.5 For every function f , 63d * df = d* d5^f + d{2HVt^ * df) — 2d{Vt^'idf). 
Proof: Let df = fiUi + /ac^a, we have *df = fiU2 - /2t^i, df = fiUJn + /2^23 and 
*df = fiuJ23 — /2t^i3- By using (jHSI) and (jHlI), we have 

h*df = {5-ifiUJ2 - hf2^i) - ^3i2df + ns[f2{auji + 6^2) - fi{buJi + CCJ2)], 
*^3df = {S3fiU2 - 5sf2UJi) - ^l3i2df + n^lfiibui - auj2) + f2{cuJi - buj2)], 
63*df - *63df = 2HVL3 *df - 2VL3Uf. (54) 

Using the operator d to act on both sides of and noticing the commutativity of d 
and ^3, we arrive at this Lemma. ^ 

Lemma 3.6 For every function f , S^d^df = d[S3{2H) * df + 2HS3 * df + 2X^3 * df — 
2HQ3*df-*d63f]. 

Proof: Similar to the proof of Lemma 13.51 we have 

h*df = 63{afi + 6/2)0^2 - ^3(6/1 + cf2)uji - Kn3 * df - 1)312^^/, 

*63df = 53(a/l + bf2)uJ2 - 63{bf, + Cf2)^l - ^Hn3 * df + Kn3 * df - 1^312 

The difference of above two equations gives 

* Ssdf -53*df = 2KVL3 *df- 2Hn3 * df. (55) 
Otherwise, It is easy to see 

*df + idf = 2H* df. (56) 

Using d63 to act on both sides of (fS^ and d to act on both sides of (j^ . considering 
the commutative relations dSs = 63d and d63 = 63d, we arrive at this Lemma. % 

If df = fiuji + f2^2, we define V/ = /iCi + /2e2, V/ = (0/1 + 6/2)61 + (6/1 + 0/2)62, 
V/ = (c/i — 6/2)61 + (a/2 — 6/1)62 and d* df = {V ■ V)fuJi A 102- It follows that, for 
function / and g, 

V/ + V/ = 2HVf (57) 
df A*dg = {Vf ■Vg)uiAu2, (58) 
rf/A*d(? = (V/- V(7)c^i Acu2, (59) 
df AUg = {S/f ■Wg)uJiAuj2. (60) 

Remark 3.1 The tensor expressions of V , V, V, V^, V ■ V, V ■ V are developed in 
[Appendix D[ 



Theorem 3.3 The second order variation of functional (jHTjl 



IS 



M 



(AH' - 2X^-^4^ - AKH-p— - 2K' 



d{2HY d{2H) OK 

dA 



M 



f)F d'^F dF 
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+ 



+ 
+ 



+ 



+ AHK 



dKd{2H) 



S + 8H' 



M 



4:{2H^ - K) 



dK 



dA 



-4- 



d£ 



dK 



d{2H)dK d{2H) 
dA 



M 



2d^£ 



d^S dS 

+ 



M 



M 



d{2HY dK 



dA 



d£ 



d{2H) 



_d{2H)dK 
d'^f ~ 2df 



V{2HQ3) ■ VQ3 



dA 



dK^ 
88 



d{2H) 



ViSH^Qs + V^fia) ■ VQ3 - V{AHQs) ■ VQ3 - 4HQ3V ■ VQ3 
- V{2Hn3) ■ Vfig - 2Hn3V ■ V^]3] dA. 



Proof: Replacing / by 1^3 in Lemma and IH.fj| and noticing that ^3 is similar to a 
constant relative to ^3, we have 

^3^ * dVL3 = [V{2Hn3) ■ Vfi3 + 2/7^3 V^fi3 - 2Vn3 ■ VO3 - 2^3^/ ■ Vn3]dA, 

63didn3 = [V(8//^fi3 + V^^]3) • Vfi3 + (8if^fi3 + V^l]3)V^fi3 

- V(4iffi3) ■ Vn3 - AHn3\/ ■ W3 - V{2Hn3) ■ Vl^g - 2Hn3V ■ Vn3]dA. 

Substituting them into (fS^ and using Lemmas l3 . 21 and we can arrive at this theorem 
through expatiatory calculations. ^ 

In particular, if dS/dK = k being a constant, ()53p is simplified to 



M 



{AH^ - 2Kf 



O^Sh d£ 



d{2Hy 



4HK 



H 



d{2H) 



+ 2K8h - 2Hp 



dA 



+ I ^^3V2n3 



M 



d{2H) 



+ 4(2if^ - K) 



d{2Hy 



H 



dA 



4d£ 



H 



mO{2H) jm 



d^£ 



H 



+ 



d£ 



H 



M 



d{2H) L 



d{2HY 

V{2HVL3) ■ VVL3 - 2VVL3 ■ VVL3 



{V^^ifdA 
dA, 



(61) 



where £h = £ — kK. 



3. 3. Shape equation of closed lipid bilayers 

Now, Let us turn to the shape equation of closed lipid bilayers. We take the free 
energy of closed lipid bilayer under the osmotic pressure as 0- Substituting £ = 
{kc/2){2H + co)^ + kK + jj, into (fS^. we obtain the shape equation Q. This equation 
is the fourth order nonlinear equation. It is not easy to obtain its special solutions. We 
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will give three typical analytical solutions as follows. Some new important results on it 
can be found in recent paper by Landolfi [22]. 

3.3.1. Constant mean curvature surface From 1956 to 1958, Alexandrov proved 
an unexpected theorem: an embedded surface (i.e. the surface does not intersect with 
itself) with constant mean curvature in must be a spherical surface [23]. Thus the 
closed bilayer with constant mean curvature must be a sphere. For a sphere with radius 
i?, we have H = —1/R and K = Substituting them into ((HI), we arrive at 

pR^ + 2/ii? - k^CQ{2 - cqR) = 0. (62) 

This equation gives the spherical radius under the osmotic pressure p. 

3.3.2. Biconcave discoid shape and \/2 torus It is instructive to find some 
axisymmetrical solutions to the shape equation 0. To do that, we denote r = 
{m cosf , M sinf , z}, ip = arctan[^^^^], and \1/ = sin-?/;. Thus (0) is transformed into 



,,2 J"^^ d'^md^i 1 fd^V p 
(\E' - 1) h \E' — - — 

du^ du^ du 2 \du J kc 



^ 2(^2 _ ^ 3^ f + (^^° + + ^ ~ 



u du^ 2u \du J \2 u kc 2m^ / du 
+ + f l)i + ll = 0. (63) 



2 kc u^ J u 2u^ 

To find the solution of ()63p that satisfies \E' = when -u = 0, we consider the 
asymptotic form of ()63|1 at u = 0: 

d^^ ^2d'^^ 1 rf^ ^ * _ ^ ^g^^ 
du^ u du^ du u^ 

Please notice that there are two misprints in our previous paper [13]. Above equation is 
the Euler differential equation and has the general solution = ai/u + a2U + a^ulnu 
with three integral constants ai = 0, a2, and a^. The asymptotic solution hints that 
= —Coulia{u/uB) might be a solution to (jUHj) which requires p = fi = 0. When 
< cqUb < e, \[' = —couln{u/uB) corresponds to the biconcave discoid shape [6,24]. 
Otherwise, when jj/kc = —2(xcq — Cq/2 and p/kc = — 2a^co, ^ = au + \/2 satisfies 
This solution corresponds to a torus with the ratio of its two radii being exactly 
V2 if a < [6,25]. 

3.4. Mechanical stability of spherical bilayers 

A spherical bilayer can be described by r = i?(sin 6* cos 0, sm6' sin 0, cos6') with R 
satisfying We have H = -1/R, K = 1/R\ V = -(l/i?)V, V ■ V = -{1/R)V^ 

and = (sin 6^) + j^,^. If we take Sh = {kc/2){2H + c,f + /i, © is 

transformed into 

-.27r 



= {2cokc/R + pR) / smOdd / #1]^ 
Jo Jo 
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+ {KcqR + 2K + pR^/2) I sm9d9 ci^fiaV^Os 



JO JO 

+ kcR^ sin Ode I d(l){V'^VL^f . (65) 

JO JO 

Expand Vt-i with the spherical harmonic functions [26]: 

oo m=l 

^3 = J2Y1 «'-^'™(^' '^)' «L = (66) 

1=0 m=-l 

If considering V^Y^m = -1(1 + ^)Yim/R^ and sin 9 d9 J^"" dcpYi^Yv^, = Smm'Sw, we 
transform into 

S^T = iR/2) J2 Wiml^W + 1) - 2]{2kjR'[l{l + 1) - coR] - p}. (67) 

l,m 

Denote that 

Pi = i2kjR^)[lil + l)-CoR] (/ = 2,3,---). (68) 

When p > pi, b^T can take negative value. Therefore, we can take the critical pressure 
as 

p, = min{pi} =p2 = i2kjR'){6 - c^R). (69) 

In this case, the spherical bilayer will be inclined to transform into the biconcave discoid 
shape. 

4. Open lipid bilayers 

In this section, we will deal with the variational problems on surface M with edge C as 
shown in FigHJ and discuss the shape equation and boundary conditions of open lipid 
bilayers with free edges. 

4-1. First order variational problems on an open surface 

In this subsection, we will discuss the first order variation of the functional 

J^= f £{2H[vlK[v])dA+ [ T{kn,kg)ds. (70) 

JM JC 

Denote = ^(2-f^[r], -^[r])(iy4, and JF^- = J^T{kn, kg)ds. 

In terms of [Appendix B[ we have uj2 = 0, ds = Ui, kn = a, kgds = Uu, Tg = b in 
curve C. Using (j2ZI)~(|SSl)! "we can arrive at 

5i^c = / d{Tn,) = 0, (71) 
Jc 



C 



'd^ 


f3T\ 


_ds^ 





,,,, dV d fdT\ d f dT 

+ 2{k^-H)kg—-rg-{—]~-(Tg — 



dkn ds \dknj ds \ dk^, 



n2ds, (72) 
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d f dV 
+ — I r. 



+ 



ds V dka 



dV 



( r —kg 1 k^ 



VLr,ds 



c \9kg 



dT 

dkry. 



kg ] ^323ds. 



In particular, (j39|) . (j49j) and (j5(J|) are still applicable. Consequently, 



d{SuJ2^ll 



M 



SuJ20,i = 0, 



c 



d{£u!iQ2 



M 



SQ2ds, 



c 



(73) 

(74) 
(75) 



M 



(V2 + 4if2 _ 2ir)-p— + (V ■ V + 2KH)^^ - 2HS 



d{2H) 



c 



+ 



62 • V 

d£ 



_d{2H) 



e2 ■ V 



dS 
OK 



dK 
d f OS 



ds \dK 



VL^ds 



dS 



Q323ds. 



(76) 



CL di2H) 

The functions ^323, ^2, and Q3 can be regarded as virtual displacements. Thus 
5jr = + 5^ + 5g)(jr^ + J^c) = gives 



(V^ + AH'^ - 2K)— + (V ■ V + 2KH)— - 2H£ = 0, 



62 ■ V 



d£ 



d{2H)_ 



d{2H) 
62 • V 



d£ 
'dK 



d 



dK 
d£ 



(77) 



c/2 f dv 



ds \ dK J ds"^ \ dkn 



dV 



dkn 



d ( dT\ d ( dV 

^^'TsK'dk^ + -(^3— 



ds V dkr, 



d£ 



d{2H) 



d£ dT dT 

kn „ _ _ ~l~ "TTZ kr, TTZ kr 



dK dK 



dk 



d^ f dT\ ^dT , , ^ 

+ K^ -k,{T- 



ds"^ \ dk^ 



d f dT 



Ta 



dk, 

d 



n 

dT_ 

dkn 



T r^, kg \ kn 
dkn 



c 



0, 



c 



(78) 
(79) 



kg + 2{kn - H)k, 



dT 
dkr,. 



dT 



£ 



0. 



c 



^ ds \dknj ds \ ^ dk, 

Among above equations, (fTTjl determines the shape of the surface M, and (f7Hj) ^ (jHn|l 
determine the position of curve C in the surface M. 



(80) 



4-.2. Shape equation and boundary conditions of open lipid bilayers 

In order to obtain the shape equation and boundary conditions of an open lipid bilayer 
with an edge C, we take £ = {kc/2){2H + cq)^ + kK + fi and T = ^kblk"^ + ^g) + 7 with 
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kf, and 7 being constants. In this case, (j77j) ^ (j8(J|) are transformed into 

K{2H + Co) (2/72 _ cqH - 2K) + kcV\2H) - 2/ii7 = 0, (81) 
kb[(fkn/ds'^ + A;„(k^/2 — rj) + Tgdkg/ds + dijgkg) / ds] 

+ to ■ V{2H) - kdTg/ds - 7A;„|^ = 0, (82) 
k^{2H + CQ) + kkn\c = Q, (83) 
kb[d'^kg/ds'^ + kg{K^ /2 — r|) — Tgdkn/ds — d(Tgkn) / ds] 

- [{kj2) {2H + cof + kK + fi + jkg] 1^ = 0, (84) 
where k"^ = k'^ + kg. 

In fact, the above four equations express the force and moment equilibrium 
equations of the surface and the edge: (jHT|) represents the force equilibrium equation 
of point in the surface M along e-^ direction; (jH^ represents the force equilibrium 
equation of point in the curve C along 63 direction; (jH^]) represents the bending moment 
equilibrium equation of point in the curve C around ei direction; ()84j) represents the 
force equilibrium equation of point in the curve C along 62 direction. 

If fcfo = 0, (jHTjl and (jH!^ remain unchanged, but (jH^ and (jH^ are simplified to 

to ■ V{2H) - kdTg/ds - 7fc„|^ = 0, (85) 
{kj2){2H + Co)' + kK + fi + -fkg\^ = 0. (86) 

4-3. Two- component lipid bilayer 

In this subsection, we study a closed bilayer consists of two domains containing different 
kinds of lipid. This problem in axisymmetrical case was theoretically discussed by 
Jiilicher and Lipowsky [27]. The shapes of two-component bilayers also were observed 
in recent experiment [28]. 

We assume that the boundary between two domains is a smooth curve and the 
bilayer is still a smooth surface. The free energy is written as 

T = p [ dV+ [ [{ki/2){2H + clf + ¥ K + ^i^]dA 

Jv J Mi 

+ / [iki'/2){2H + c'o')'' + k''K + fi'']dA + ^ [ ds. (87) 

Jmii Jc 

In terms of the discussions on closed bilayers in section El and above discussions in 
this section, we can promptly write the shape equations of the two-component bilayer 
without any symmetrical assumption as 

p - 2fi'H + ki{2H + Co){2H^ - d^H - 2K) + ky\2H) = 0, (88) 

where the superscripts i = I and II represent the two lipid domains, respectively. And 
the boundary conditions are as follows: 

fcfes ■ V(2if) - VdTg/ds + kl}e2 ■ V{2H) - k^^drg/ds - 7A;„|^ = 0, (89) 
ki{2H + c^) + k'kn - [ki\2H + c'J) + V'knWc = 0' (90) 
{ki/2)i2H + cif + ¥K + ^l' 
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- [{ki'/2){2H + ci'f + k"K + + 7A;^|^ = 0. (91) 

In above equations, the positive direction of curve C is set to along ei of the hpid 
domain consisting of component I. Furthermore, ()88|) ^ fj91|) are also applied to describe 
the closed bilayer with more than two domains. But the boundary conditions are not 
applied to the bilayer with a sharp angle across the boundary between domains. 



5. Cell membranes with cross-linking structures 

Cell membranes contain cross-linking protein structures. As is well known, rubber also 
consists of cross-linking polymer structures [29]. In this section, we first drive the free 
energy of cell membrane with cross-linking protein structure by analogy with the rubber 
elasticity. Secondly, we derive the shape equation and in-plane stain equations of cell 
membrane by taking the first order variation of the free energy. Lastly, we discuss the 
mechanical stability of spherical cell membrane. 



5.1. The free energy of cell membrane 

Above all, we discuss the free energy change of a Gaussian chain in a small strain field 

(92) 



-yy 



\ 







\ 





with Ezz = —{^xx + £xx) expressed in an orthogonal coordinate system Oxyz. 

Assume that one end of the chain is fixed at origin O while another is denoted by 
Riv before undergoing the strains, where N is the number of the segments of the chain. 
The partition function of the chain can be calculated by path integrals [30]: 



Ro 



D[Rn] exp 



2L2 



N 



dn 



gR„ 

dn 



aexp 



3(R 



■N 



Ro)' 



2iVL2 



where a is a constant, and L is the segment length. After undergoing the strains, the 
partition function is changed to 



cxexp 



3(R 



Rn 



2iVL2 



Considering the relation (Rat — Ro)^ 
function of end-to-end distance P{Rn — 
change as a result of the strains: 



= [(1 + e) ■ (Rtv — Ro)]' and the distribution 
Ro) = J dRj^z ' obtain the free energy 



/. = - ksT J dRr,{\nZ, - lnZ)P(R^ - Rq) 

= ksTKexx + £yy) ~ i^xx^yy ~ ^xy)]- 



'xy/\- (93) 

For the cell membrane with the cross-linking structure, we assume that: (i) The 
membrane is a smooth surface and junction points between protein chains are confined in 
the vicinity of the surface and freely depart from it within the range of ±/i/2, where h is 
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the thickness of the membrane, (ii) There is no change of volume occupied by the cross- 
hnking structure on deformation. On average, there are A/" protein chains per volume, 
(iii) The protein chain can be regarded as a Gaussian chain with the mean end-to-end 
distance much smaller than the dimension of the cell membrane, (iv) The junction 
points move on deformation as if they were embedded in an elastic continuum (Affine 
deformation assumption), (v) The free energy of the cell membrane on deformation is 
the sum of the free energies of closed lipid bilayer and the cross-linking structure. The 
free energy of cross-linking structure is the sum of the free energies of individual protein 
chains. 

If the cell membrane undergoing the small in-plane deformation I | , 

\ £12 £22 J 

where "1" and "2" represent two orthogonal directions of the membrane surface, we 
can obtain the free energy of a protein chain on deformation with the similar form of 
fl^ under above assumptions (i)~(iv). Using above assumption (v), we have the free 
energy of the cell membrane under the osmotic pressure p: 



T= [ {Sd + SH)dA + p ! dV, 

J M Jv 



where Sh = {kc/2){2H + cq)^ + /i and Sd 

2 J = £11 + 622, Q = £\\£22 ~ £1 



{ka/2)[{2Jf 



(94) 

Q] with kd = 2XhkBT, 



-2 
-12- 



Remark 5.1 We do not write the term of kK in fl94|) because j^^kKdA is a constant 
for closed surface (see [Appendix C ). 



5.2. Strain analysis 

If a point fq in a surface undergoing a displacement u to arrive at point r, we have 
d\i = dr — dro and Sidu = 6idr {i = 1, 2, 3). 

If denote dr = UiSi + 1^262 and du = Uicui + U2i-o'2 with |Ui| ^ 1, IU2I ^ 1, we 
can define the strains [31]: 



£11 



£^22 



du ■ ei 
\dro\ 

du ■ 62 
\dro\ 



UJ2=0 



wi=0 



eu-- 



du ■ e2 



Ui ■ ei, 

U2 • 62, 

du - Bi 

\dro\ 



(95) 
(96) 



-(Ui-e2 + U2-ei). (97) 



(^2=0 \ l""-U| / a;i=0 

Using 6idu = 6idr and the definitions of strains (j95|) ^ (j97j) . we can obtain the 
variational relations: 

diEiiUJi A CJ2 = (1 - £^ii)5i^^i A CJ2 - U2 ■ ei6iUJ2 A UJ2 

+ fiii2Ui • e2UJi ALO2 + fijisUi • e^LOi A 102, 

SiEuUJi A CJ2 = ^[(1 - eu)uji A 6iUJi + (1 - e22)SiUJ2 A u;2 - U2 ■ eiUi A SiUJ2 
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- Ui ■ e25iUJi Auj2 + fij2i(^ii - £^22)^^1 A UJ2 
+ fii23Ui ■ e^ui Au2 + ^ii3U2 ■ escul A U2], 
5ie22^i A 0^2 = (1 — £^22)^^1 A 5iUJ2 — Ui ■ e2UJi A 5iUJi 

+ fii2iU2 ■ eicji A u;2 + fij23U2 • escui A uj2- 
The leading terms of above relations are: 

SiEiiUJi A Co'2 = 5iUJi A UJ2, 

6iei2UJi Auj2 = -[uJi A 5iUJi + 5iUJ2 A ^^2], 



(5ie22t^l A CJ2 = CJi A 5iUJ2. 



Thus, 



5i{2J)uJi Auj2 = 5i{eii + £22)^^^! A ti;2 = 5iUJi Auj2 + uJi A 5iUJ2, 
6iQuJi Auj2 = 6i{eu£22 - £^12)^1 A ^2 

= (eiicui + 6:12^^2) A (5ju;2 - (^12^^! + £22^2) A SiUJi. 
Considering equations (I2Z|)~(ISS1), (jl01|) and ()102p . we have 
(5i(2J)c<Ji A co'2 = (i(i7iCt;2), 

^iQcui A CJ2 = {£iiduj2 - ei2duJi)Q,i - {euloi + 622^^2) A d^li, 
62{2J)uJi A 0^2 = — d{Q2^^i), 

62QUJ1 Auj2 = (eiicc'i + ei2C^2) A d^l2 + ^2{£i2duj2 - e22duJi); 
63{2J)uJi Auj2= - 2HnsdA, 

SsQuJi Auj2 = [-2H{2J) + aeii + 2hei2 + ce22]^3dA. 



(98) 
(99) 
(100) 

(101) 

(102) 

(103) 
(104) 

(105) 
(106) 

(107) 
(108) 



5.3. Shape equation and in-plane strain equations of cell membranes 

To obtain the shape equation and in-plane strain equations of cell membranes, we 
must take the first order variation of the functional (j94j) . Denote Td = JM^^dA and 
^cp = Jm ^udA + pjy dV. 

From p03|) ^ ()108|) . we can calculate that: 



M 



■^^6i{2J)dA+-^SiQdA + Fdi2J,Q)S,dA 



d 



M 



' M 



BF 1 C BF 

—^ + Fd{2J,Q) Au;2fii+/ ^{€uduj2 - ei2duji)ni 

dFd 



(^12'^l + £^22'^2, 



dQ 



(109) 



52^d 



M 



BF BF 
— 4^52(2J)rfA + ^52QdA + Fd{2J, Q)62dA 
B[2J) BQ 



d 



M 



dFd 
_B{2J) 



+ £di2J,Q) 



A UJ1Q2 



dFd 



Qod 



M 



dFd 
dQ 



[{snUJi + 612^^2) 



{ei2duj2 — e22duji)n2 
(110) 
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M 



M 



dQ 



[aeii + 2bei2 + 0^22] ^sc?^ 



- / 2H 



VL^dA. 



[Ill) 



Otherwise, section El tells us: 



cp 



0, 



~p 



M 



(V^ + m"^ - 2K) 



dS 



H 



d{2H) 



2H£h + p 



Q'.dA. 



Therefore, 6iJ-' = 6iJ-'d + SiJ-'cp = gives: 



d 



_d{2J) 



- d 

dSd 
d{2J) 

+ d 



(£121^1 + 822^2 



dQ 

d£d 
dQ 



0, 



(112) 



d£d 
dQ 



d£ 

^ ^1 + ^i^i2duj2 - e22duji] 
oQ 



{eiiui + ei2uj2) 



(V^ + AH^ - 2K) 
dS. 



d£ 



H 



d{2H) 



2H 



0, 



(113) 



dSd 



d{2J) 



£d + (2 J) 



d£d 
dQ 



+ P+ "^["^11 + 2^^i2 + ce22] = 0. 



(114) 



Substituting £h = ^{2H + cq)^ + /i and £d = ^[{2JY — Q] into above three equations. 



we obtain: 



kd[-d{2J) Auj2- ^{eiiduj2 - e^dui) + ^d{ei2UJi + 622^^2) 



0, 



kd[d{2J) A c^i - ^{ei2duj2 - £22^^^!) - ^d^euUJi + 612^^2)] = 0, 

p - 2H{fi + kdJ) + kc{2H + co)(2if2 - coH - 2K) + kcV\2H) 
kd 

- y (ae^ii + 26£:i2 + 0^22) = 0. 



(115) 
(116) 

(117) 



flllSj) and piCij) are called in-plane strain equations of the cell membrane, while flll7j) 
is the shape equation. 

Remark 5.2 The higher order terms of Sij (i,j = 1,2) are neglected in above three 
equations. 

Obviously, if fc^ = 0, then pi5j] and pifij) are two identities. Moreover pi7j] 
degenerates into shape equation (jS)) of closed lipid bilayers in this case. Otherwise, for 
small strain, ()117j) is very close to (jH)), which may suggest that cross-linking structures 
have small effects on the shape of lipid bilayers. 



Geometric theory on the elasticity of bio-membranes 



20 



It is not hard to verify that £12 = 0, en = 622 = £ satisfy (jll5|) and if e being 
a constant. In this case, the sphere with radius R is the solution of ()117p if it satisfies 



5.4- Mechanical stabilities of spherical cell membranes 

To discuss the stabilities of spherical cell membranes, we must discuss the second order 
variations of the functional JF. In mathematical point of view presented in section |21 we 
must calculate S^SjJ-" {i, j = 1, 2, 3). But in physical and symmetric point of view, we just 
need to calculate because we can expect that the perturbations along the normal 
are primary to the instabilities of spherical membranes under the osmotic pressure which 
is perpendicular to the sphere surfaces. 



If we taking Sh = ^{2H + cof + /x and Sd = ^[(2-^)^ - Q], the leading term of 
(UnD is 



pR^ + (2/i + 3kd£)R + hcoicoR - 2) = 0. 






Using Lemmas EH E21 Eqs.dHHl^ (11001), we can obtain 



^l^d = - ^ / [("^^a^ii + 253^12 + 53C£22) + 63{2H){2J)] n^dA 



Jm 




(120) 



for the spherical cell membrane with radius R and strain e. 
Otherwise, (jF)T|l suggests that 




Therefore 



51:f 



^3^d + S^^cp 




+ / nsV'^ns{2k^co/R + 2k J R^ - {SkdS + k^cl + 2/i) /2}dA 



Jm 




(122) 
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If considering (jll8p and expanding ^3 as (jUUI) . we have 
SlJ^= j nl{3ka/R^ + i2Kc^/R^) +p/R}dA 

+ [ n3V^n3{k,co/R + 2kjR^ +pR/2}dA+ I k^iV'^Vt^fdA 
Jm Jm 

= Yl +[/(/ + 1) - 2][l{l + l)kjR^ - k,co/R - pR/2]}. 



l,m 

|2 



The zero point of the coefficient of |aim| above expression is 

n - 2kM + ^)-CoR] . . . . 

^' - [/(/ + 1) - 2]R + R^ (^-2,3,---). (123) 

Obviously, on the one hand, if A;;^ = 0, (jl23p is degenerated into (jUH)) with I > 2. On 
the other hand, if fc^ > 0, we must take the minimum of ()123p to obtain the critical 
pressure. 

If let ^ = /(/ + 1) > 6, we have 

dp Qkfi/R 2kc 

di = ~J^^^R^' ^^^^^ 

dp/d^ = and ^ > 6 imply ^ = 2 + R^J'ikd/kc which is valid only if 3kdR^ > IQkc- 
Therefore, the critical pressure is: 

r 3fcrf I 2fcc(6-coR) 2fce(10-coR) CjU jd2 ^ ^ ci, \ 

^^"""^^'^"1 ^ + t(2-coi?) (3^./?^ > 16fc.). ^''^^ 

Eg . (11271) includes the classical result for stability of elastic shell. The critical 
pressure for classical spherical shell is pc oc Yh'^/R^ [32,33], where Y is the Young's 
modulus of the shell. If taking cq = 0, fc^ oc Yh^ kc oc Yh'^ and R^ h, our result ()127|) 
also gives Pc oc Yh'^/R^. As far as we know, this is the first time to obtain the critical 
pressure for spherical shell through the second order variation of free energy without 
any assumption to the shape of its losing the stability (cf. Ref. [32,33]). 

Otherwise, if we take the typical parameters of cell membranes as kc ~ 20kBT [4,5], 
kd ~ 2A^N/m [34], h ~ 4nm, R ~ 1/im, cqR ~ 1, we obtain ~ 4 Pa from (|127p . 
which is much larger than p^ ~ 0.2Pa without considering k^ induced by the cross- 
linking structures. Therefore, cross-linking structures greatly enhance the mechanical 
stabilities of cell membranes. 



6. Conclusion 



In above discussion, we deal with variational problems on closed and open surfaces by 
using exterior differential forms. We obtain the shape equation of closed lipid bilayers. 
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the shape equation and boundary conditions of open hpid bilayers and two- component 
hpid bilayers, and the shape equation and in-plane stain equations of cell membranes 
with cross-linking protein structures. Furthermore, we discuss the mechanical stabilities 
of spherical lipid bilayers and cell membranes. 
Some new results are obtained as follows: 

(i) The fundamental variational equations in a surface: Eqs. (fTfj) ~ (jHBj) . 

(ii) The general expressions of the second order variation of the free energy for 
closed lipid bilayers: theorem |H.HI and Eq. (pTT|) . 

(iii) The general shape equation and boundary conditions of open lipid bilayers and 
two-component lipid bilayers: Eqs. (fTTjl ~(jHni) and (jHH|) ^ (jnT|l . 

(iv) The free energy (jMjl . shape equation and strain equations pi5|) ^ pi7|) of the 
cell membranes with cross-linking protein structures. 

(v) The critical pressure ()127j) of losing stabilities for spherical cell membranes. It 
includes the critical pressures not only for closed lipid bilayers, but also for the classic 
solid shells. Otherwise it suggests that cross-linking protein structures can enhance the 
stabilities of cell membranes. 

In the future, we will devote ourselves to applying above results to explain the 
shapes of open lipid bilayers found by Saitoh et al, and predict new shapes of multi- 
component lipid bilayers and cell membranes. Moreover, We will discuss whether and 
how the in-plane modes affect the instability of cell membranes although we believe they 
have no qualitative effect on our results in section 5.4. 
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Appendix A. Exterior differential forms and Stokes' theorem 

A manifold can be roughly regarded as a multi-dimensional surface. In the neighborhood 
of every point, we can construct the local coordinates {u^, u^, ■ ■ ■ , u"*), where m is the 
dimension of the surface. In this paper we just consider smooth, orientable manifolds 
and smooth functions. 

We call the function /(m^, m^, ■ ■ ■ , m"*) 0-form and aj(M^, u^, ■ ■ ■ , 1-form, 
where Einstein summation rule is used and it is also used in the following contents. The 
r-form (r < m) is defined as ai^i^.-i^du^^ A dw*^ A ■ ■ ■ A du''' , where the exterior production 
"A" satisfies du' A du^ = —du^ A du'. Denote A'' = {all r-forms}, (r = 0, 1, 2, ■ ■ ■ , m). 
Definition A linear operator d : A'' — > A''^^ is called the exterior differential operator 
if it satisfies: 

(i) For function f{u^, u^, ■ ■ ■ , m"*), df = §^du'' is an ordinary differential; 
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(ii) dd = 0; 

(iii) WuJi e A'' and Vco'2 e A'', d{uJi A LJ2) = dui Alj2 + {-IYlji A duj2. 

Stokes theorem If u is an (m — l)-form with compact support set on M, and P is a 
domain with boundary &D in M, then 

duj= 00. (A.l) 

V JdV 

Appendix B. Curves in a surface 

If a curve passes through P in the surface, we construct a Frenet frame {T, N, B} such 
that T, N and B are the tangent, normal and binormal vectors of the curve, respectively. 
Denote 6 the angle between ei and T. Set M = es x T. Thus we have 

ei cos 6* + e2 sin 6*, 
= — ei sin 6 + 62 cos 6. 

It is not hard to calculate 

dT = {de + M + eg (c^ig cos 9 + uj2z sin 9) . 

Frenet Formulas tell us dT/ds = kN. Therefore, we have the geodesic curvature, 
the geodesic torsion, and the normal curvature of the curve: 

kg = kN-M = (dT/ds) ■M= {d9 + UO12)/ ds, (B.l) 
Tg = —(de-^/ds) ■ M = [{bui + CUJ2) cos 9 — {aui + bu2) {sin 9)] /ds, 
kn = II /I = (at^i + 2buJiU!2 + cuD/iojl + uD- 
If the curve along ei such that ^ = 0, we have ds = uji,uj2 = Q and 

kg = uju/uji, Tg = b, and fc„ = a. (B.2) 

Appendix C. Gauss-Bonnet formula 

Using ©, (Uni) and (HSI), we have 

dLJi2 = —KuJi A UJ2. (C.l) 

This formula was called Theorem Egregium by Gauss. From Theorem Egregium and 
()B.1|) . we can derive Gauss-Bonnet formula: 

! KdA+ I kgds = 2nx{M), (C.2) 
Jm Jc 

where x{^) is the characteristic number of smooth surface M with smooth edge C. 
x(M) = 1 for a simple surface with an edge. For a closed surface, we have 

I KdA = 2-Kx{M). (C.3) 

J M 
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Appendix D. The tensor expressions of V, V, V, V^, V • V, and V • V 

At every point r in the surface, we can take local coordinates {u^,u'^) where the first 
and the second fundamental form are denoted by / = gijdu^du^ and // = Lijdu^du^ , 
respectively. Let (g^^) = {gij)~^, (U^) = {Lij)~^ and rj = dr/du\ thus we have 

d 



V = Yi{2Hg'^ - KV^) 



d 



V- V 



V- V 



_]_d_ 

1 



d 



\f99 



d 



d 



Vgdu' 

J__d_ 

y/gdu' 



Vgi^Hg^^ 



d 



d_ 



As an example, We will prove the last one of above expressions. 
Proof: If taking the orthogonal local coordinates, we have / = gii{du^Y + 922{du'^Y — 
+ UJ2, which implies uji = ^gudu^ and 002 = ^/go-idu^. For function /, on the one 
hand, we have df{u^,u^) = fiUi + f2U!2 = fiJlhidu^ + f^Jlfndu^, on the other hand. 



we have df = S^du' + ^du\ Therefore, h = f2 = 



The second fundamental form II — acuf + 2huJiU}2 



322 dv? ■ 



a = Lii/gu, b = Lvilsfg, c 
and 

L22 



ll-t^22 



12 



L 

^22 



-^12 

n2 



L11L22 — Lu 
Ln 



L11L22 — L 



1 

-^22/5'22- Thus K ^ ac - b'^ 

>L22 = gKL'^- 

Lu = -gKL'^; 
Ln = gKL^^. 



Lijdu^du' implies 
(L11L22 — L\^/g^ 



12 



Moreover, we have 



d*df 



/2'^13 + /l'^23 = -f2{aUJi 



1 



L12— - i-ii— I du + 



6CJ2) 
1 



d 
+ 



d 



Therefore, V • V/ 



_ j_ 

u)i/\W2 a/5 



du 

V9K 



^ dv? ) 

11 df 



V9 



f2{buji - 
dl 



CUJ2) 



L 



22 



-L 



12 



dl 
du^ 



du^; 



+ L 



12 



df 



du' 

rl2 df 



du^ 



du 



+ 1^- 



22 



dl 
du^ 



du^ A du^. 
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